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Common origin of the high energy astronomical 
gamma rays, neutrinos and cosmic ray positrons? 

Shlomo Dado^ and Arnon Dar^ 

ABSTRACT 

We show that the observed fluxes, spectra and sky distributions of the high 
energy astronomical neutrinos, gamma rays and cosmic ray positrons satisfy the 
simple relations expected from their common production in hadronic collisions 
in/near source of high energy cosmic rays with diffuse matter. 

Subject headings: (ISM): gamma rays, neutrinos, dark matter 


1. Introduction 

Nearly half a century ago it has already been realized that the decay of mesons produced 
in the interactions of high energy cosmic rays (CRs) with matter and radiation in/near the 
cosmic ray sources, in the interstellar medium (ISM), and in the intergalactic medium (IGM), 
can produce detectable fluxes of high energy gamma rays ( 7 ’s), neutrinos (z^’s) and positrons 
(e“''’s), which carry unique information on the ultra relativistic universe (e.g., Zatsepin & 
Kuzmin 1966; Berezinsky & Smirnov 1975; Margolis et ah 1978; Eichler 1978; Eichler & 
Schramm 1978; Stecker 1979). In particular, interaction of cosmic rays with diffuse matter 
and/or radiation in/near the cosmic ray sources, in the ISM of Galaxies and the IGM has 
long been recognized as possible point and diffuse sources of very high energy Galactic 
and extragalactic 7 -rays and neutrinos (see, e.g., Dar 1983,1985,1991; Halzen et al. 1990; 
Berezinsky 1991; Stecker et al. 1991; Mannheim 1993; Berezinsky et al. 1993,1994; Domokos 
et al. 1993; Dar & Shaviv 1996; Halzen 1996; Dar & Laor 1997). 

The successful detection of neutrinos from supernova 1987A with the deep underground 
Kamiokande and IMB water Gherenkov telescopes, of Galactic and extragalactic point 
sources of very high energy 7 -rays with ground based atmospheric and water Gherenkov 
telescopes (see, e.g., Weekes et al. 1989; Punch et al. 1992), and of point sources and 
diffuse Galactic and extragalactic backgrounds of high energy 7 -rays with the Gompton 
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Gamma Ray Observatory (Hunter et al. 1997; Sreekumar et al. 1998), have led to the de¬ 
velopment of new generations of much larger and more sensitive telescopes of high energy 
7 -rays, neutrinos and cosmic rays, which have already made several breakthrough obser¬ 
vations. Such breakthroughs include the discovery of a diffuse high energy neutrino back¬ 
ground radiation (NBR) with the megaton IceCube detector under the south pole (Aartsen 
et al. 2013,2014,2015), the accurate measurement of the diffuse high energy gamma-ray back¬ 
ground radiation (GBR) with the Large Area Telescope (LAT) aboard the Fermi satellite 
(Ackermann et al. 2012,2015), and the precise measurements of the high energy local back¬ 
ground of cosmic ray positrons with the Payload for Antimatter Matter Exploration and 
Light-nuclei Astrophysics (PAMELA) satellite (Adriani et al. 2013) and the Alpha Magnetic 
Spectrometer (AMS) aboard the International Space Station (Aguilar et al. 2014). 

Possible sources and the detailed properties of the observed astronomical neutrinos, 
7 -rays and positrons were discussed in the discovery papers and in a large number of pub¬ 
lications, which followed the observations (for a recent review see, e.g., Anchordoqui et 
al. 2014a). But, most of the published theoretical estimates of the fluxees of the high en¬ 
ergy neutrinos, 7 -rays and positrons from point sources and diffuse sources, both before and 
after these recent breakthrough observations, were fraught with uncertainties concerning the 
production mechanisms of the high energy neutrinos, 7 -rays and positrons and the proper¬ 
ties of their sources, which have left the origin and observed properties of these radiations 
essentially nnsolved puzzles. 

In view of snch dihculties, several anthors have tried alternative approaches towards re¬ 
solving these puzzles, such as deriving npper bonnds on the observed fluxes of the high energy 
astronomical neutrinos, gamma-rays and positrons, or general relations between them, which 
may indicate their origin. E.g., a highly cited ’’model-independent npper bonnd” on the en¬ 
ergy flux of extragalactic NBR above 100 TeV was derived by Waxman & Bahcall (1999). It 
was based on the assnmption that the NBR is prodnced by the decay of vr^'^’s prodnced by the 
extragalactic cosmic ray protons in proton-photon collisions in the extragalctic GR sonrces. 
However, as shown in Appendix Al, the ’’model independent npper bonnd” of Waxman & 
Bahcall (1999), actnally, is a model dependent estimate based on an incorrect extrapolation 
of the observed flnx of extragalactic GRs to lower energies and on incomplete conseqnences 
of their arbitrary choices. When properly corrected, this ’’npper bonnd” increases by nearly 
two orders of magnitude! (see Appendix Al) and becomes useless. 

Other authors have tried to relate the NBR discovered with IceGube and the high energy 
gamma-ray backgronnd radiation (GBR) measnred with Fermi-LAT, assnming that both 
were produced in hadronic interactions of high energy cosmic rays (e.g., Mnrase et al. 2013; 
Anchordoqni et al. 2014a and references therein). However, most of these attempts were 
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based on the uncertain assumptions that the NBR is mostly extragalactic in origin and that 
the very high energy 7 -rays from blazars and other extragalactic point sources are leptonic 
in origin. Consequently, the NBR was related only to the diffuse intergalactic gamma-ray 
background (IGRB), ignoring the possibility of a much larger Galactic contribution and an 
extragalactic contribution (e.g. from blazars): 

Multiwavelength observations of blazars have shown that their spectral energy density 
usually exhibits two bumps, one at low energy (infrared to X-rays) and one at high energy 
(see, e.g., Abdo et al. 2010 and references therein). The origin of the high energy bump near 
TeV is still under debate (see, e.g., Cerruti et al. 2015 and references therein). Moreover, 
current Cherenkov telescopes have identified a population of ultra-high-frequency peaked BL 
Lac objects (UHBLs), which exhibit exceptionally hard TeV spectra. This hard emission, as 
well as the high energy emission in gamma-ray bursts (GRBs) challenge the synchrotron-self- 
Compton (SSC) model (e.g., Konigl 1981) of very high energy 7 -ray emission from blazar 
jets (Maraschi et al. 1992; Bloom & Marscher 1996) and from GRB jets (e.g.. Dado & 
Dar 2005,2009). In particular, the standard one zone synchrotron-self-Compton (SSC) in¬ 
terpretation of the double peaked blazer spectra (see e.g., Abdo et al. 2010 and references 
therein) implies that the SSC peak frequency is inversely proportional to the synchrotron 
peak frequency because of the Klein Nishina suppression of the high energy inverse Compton 
scattering cross section. Such a correlation does not seem to be present in the observational 
data. In fact, the Klein-Nishina suppression of inverse Compton scattering of background 
photons by high energy electrons and the suppression of the population of high energy elec¬ 
trons by synchrotron radiation and inverse Compton scattering from background photons 
suggest that 7 r° —)■ 2 7 decay following inclusive production and decay of mesons in jets en¬ 
counters with diffuse matter may become the dominant production mechanism of very high 
energy photons emitted by blazars (Dar & Shaviv 1996; Dar & Laor 1997) and GRBs (Dado 
& Dar 2005,2009). 

If the main origin of the very high energy astronomical neutrinos and 7 -rays and CR 
positrons observed near Earth is the decay of mesons produced in high energy CR collisions 
with diffuse matter in/near the CR sources, and/or in the ISM, and/or in the IGM, then 
under very general assumptions their fluxes, spectra and sky distributions are simply related. 
In this paper, we reconsider such relations (e.g.. Dado & Dar 2014) and confront them with 
observations rather than prejudices. Indeed, we show that the high energy NBR observed 
with IceCube and the GBR observed with Fermi-LAT (and the high energy Cherenkov 
telescopes) satisfy the simple relations expected from hadronic production of mesons by 
high energy cosmic rays in collisions with diffuse matter in/near source (rather than with 
diffuse radiation - see Appendix A2). They imply a non isotropic NBR whose sky distribution 
is similar to that of the Fermi-LAT GBR. We also show that the flux of the high energy CR 
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positrons observed near Earth with PAMELA and AMS02 is that expected from secondary 
prodnction of mesons in the local ISM by the flnx of cosmic ray nncleons (protons and 
nncleons bonnd in atomic nnclei) observed near Earth. 


2. CR production of secondaries 

Hadronic collisions of high energy nncleons of energy Ep with diffnse matter prodnce 
7 -rays, z/’s and e^’s mainly throngh tt and K decays. If the lab frame energy E of these 
secondary particles is expressed as xEp, the distribution of x is independent of Ep to a good 
approximation (Feynman 1969). Consequently, a flux <hp oc E"^ of CRs in a diffuse matter 
produces through hadronic collisions secondary y-rays, z/’s, and e^’s with a flux per unit 
volume (e.g., Dar 1983) 

<I>i(E) oc ain(Ep) Us c <l>p(Ep)/xi (1) 

where i = 7 , z/, or e^, the mean value x, depends only on the distribution of x in the inclusive 
production pp —)■ iX but not on E, Ep = E/x, is the mean energy of cosmic ray protons that 
produce particles i with energy E, and ain ~ 30 x (Ep/GeV)°'°^® mb is the pp total inelastic 
cross section. 

Eq. (1) implies that CR collisions with diffuse matter produce secondary fluxes of 7 - 
rays, z/’s and e^’s with the same power-law spectrum ~ g-z^+o-oss different normalization. 
These fluxes are later modihed by propagation effects: attenuation of high energy 7 - rays 
mainly by Compton scattering from free electrons and by pair production on background 
photons, oscillations of neutrinos in space that spread the neutrino flux over the three neu¬ 
trino flavors, diffusion propagation of high energy e^’s in magnetic helds, and energy loss of 
e^’s by synchrotron emission and inverse Compton scattering of background photons. 


3. The GBR-NBR connection 

We shall assume that the main kinds of sources of high energy cosmic rays (supernova 
explosions, merger of compact stars, and mass accretion onto compact stars, stellar mass 
black holes and massive black holes) are common to our Galaxy and external galaxies, that 
Fermi acceleration (Fermi 1949) and galactic conhnement produce a power-law spectrum 
<hp ~ E“^ of galactic CRs, and that the main source of background 7 -rays and neutrinos 
above TeV is hadronic collisions of these high energy CRs with diffuse matter inside or near 
the CR sources, in the ISM of their host galaxies, and in the IGM, which produce mesons 
(mostly pious and K mesons) that decay to high energy 7 -rays, z/’s and e^’s. The assumption 
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of hadronic rather than leptonic production is based on the facts that inverse Compton 
scattering of background photons by high energy electrons is suppressed by the Klein Nishina 
cross section and that the population of high energy electrons is suppressed by synchrotron 
radiation and inverse Compton scattering from background photons. (Photoproduction has 
an effective threshold too high to explain the high energy GBR observed below TeV). The 
hadronic production hypothesis, mostly overlooked, will be confronted here with observations 
rather than with prejudices. 

Consider hrst hadronic production cf high energy 7 -rays and neutrinos by cosmic ray 
collisions with diffuse matter inside or near the CR sources. 

The local flux of high energy cosmic ray nucleons (free protons and nucleons bound in 
atomic nuclei that hereafter will be denoted by p) between several GeV and the cosmic ray 
knee energy Eknee ~ 1 PeV/nucleon is well described by 

$p(E) ^ C(E/GeV)-^ fu (2) 

where C ks 1.8, (3 ^ 2.70 (e.g., Olive et al. 2014), and fu = l/(GeV cm^ s sr) is the flux unit. 
Between the CR ’knee’ and CR ankle at E~ 3 EeV, C ~ 114, and jS ~ 3.0 (Apel et al. 2013). 

We will ignore the small differences between CR protons and nucleons bound in CR nuclei 
(A,Z) at the same energy per nucleon, because such nucleons contribute only a few percents 
to the total CR p flux and most of them are bound in light nuclei whose inelastic cross 
section per nucleon is ~ (Tpa/A Upp, i.e., roughly the same as that of free protons. 

CR protons escape the Galaxy by diffusion through its turbulent magnetic helds. For 
a Kolmogorov spectrum (Kolmogorov 1941) of random magnetic helds, their escape time 
satishes resc(E) oc In a steady state, the supply rate of high energy CR protons by 

the Galactic CR sources (s) per unit volume is equal to their escape rate from the Galaxy. 
Hence, the source spectrum of CR protons satishes <hs oc <hp/resc oc E“^" where /Sg = /9 — 1/3. 
Roughly, /9s = /9 — 1/3;^ 2.37 for E < Eknee (p), which is consistent with Fermi acceleration 
modihed by escape by dihusion, and /Sg ~ 2.67 for Eknee(p) < E < Eankie(p)- 

We shall assume that above 100 GeV the contribution of electron bremsstrahlung and 
inverse Compton scattering of background photons to the GBR, that decreases with in¬ 
creasing energy, becomes relatively small compared to the contribution from Galactic and 
extragalactic 7 -ray production by 7 r° 27 decays. Neutrinos are produced mainly by tt and 

K decays. 

As long as the neutrinos and gamma-rays are produced mainly through tt decay by CRs 
with energy below the CR knee whose in/near source spectral index is /9j ~ 2.7 — 1/3 — 0.06 ~ 2.31, 
their huxes satisfy 


<h,y(E) (m^+/2mjro)^'^^ ‘^ 7 (E) ~ 0.39<h.y(E) 


( 3 ) 
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per V flavor. At higher energies neutrinos are produced mainly through —)■ z/ and 

K+ —)■ decays followed by —)■ e^ 3z/ decay, and in —)■ 7r’''7r°; 27r+, tt”; 27r° 7r+; 7r° /i’*' 

and 

K° —)■ 7 r’'' 7 r“; 7 r° 7 r°; 7 r+ 7 r“ 7 r°; decays followed by the and /i^ decays. The 7 -rays are 

produced mainly by 7 r° —)■ 27 decays. Well above the CR knee where (3^ ~3 — 1/3 — 0.06 ~ 2.61, 
inside/near source hadronic meson production by cosmic rays yields 

$^(E) ^ 0.52$..^(E). (4) 

per z/ flavor. Eqs. (3),(4), cannot be tested directly with current data because they neglect 
the attenuation of high energy 7 -rays, and because the NBR was measured at E^, > 30 TeV, 
while the published Fermi-LAT data on the GBR (Ackermann et al. 2012) and on the 
extragalactic gamma background (Ackermann et al. 2015) are limited to E< 1.2 TeV (see 
Figs. 1,2) (Other recently reported measurements of the high energy GBR with Gherenkov 
telescopes such as those with H.E.S.S around 15 TeV (Abramowski et al. 2014) and in 
the ARGO-YBJ experiment around 1 TeV, (Ma 2011) were limited to the Galactic plane 
|b| < 2 °.) 

The attenuation of Galactic 7 -rays of energy below TeV is negligible, while the observed 
extragalactic 7 -ray background (EGB) is strongly absorbed above 100 GeV. The full sky 
GBR measured with Fermi-LAT in the energy range 100 GeV - 2 TeV can be corrected 
for the attenuation of the EGB and used in Eq.(3) to predict in this energy range, and 
then it can be extrapolated to the energy range of the NBR detected with IceGube using 
4>u{^) cc (Jin((l-l-z) 20 E) <l>s((l-|-z) 20 E), which follows from Eq. (1) for extragalactic GR 
sources. In this relation, l-|-z^ 2.5 ± 0.5 is the redshift at the peak of star formation rate 
(i.e., of supernova explosions and GRBs) and of the evolution function of the emission by 
BL Lac objects - presumably the main extragalactic sources of high energy GRs. 

Eq. (1) predicts that the non attenuated EGB behaves like well below the GR 

knee. Indeed, the EGB measured with Fermi-LAT below 820 GeV was best ht (Ackermann 
et al. 2015, Model G) with an exponential cutoff power-law, 

$EGB ~ (6.42 ± 0.40) X 10"^ (E/GeV)-2-^°^°-°2 fu (5) 

where Ec ~ 366 ± 100 GeV (y^/df = 6.9/23). This best ht is shown in Fig. 1. Presumably, 
the power-law represents the unattenuated EGB at energies well below the ’’knee”, produced 
by high energy cosmic rays in external galaxies. Eqs. (3) and (5) then yield (piy ~ 1.03 x 10“^^ 
fu per u havor at E=100 GeV whose extrapolation to E> 50 TeV yields an isotropic extra¬ 
galactic [EG] neutrino hux per z/ havor, 

r 1 -0.61±.05 

E2<I>,[EG] ^ (0.85±.30) X 10-® GeV^fu. ( 6 ) 
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Similarly, the Galactic [MW] contribution (^gbr — ‘hEGe) ~ 2.49 x 10“^ GeV^ fu 

to the GBR at E=TeV and Eq. (3) can be used to estimate <hi,[MW], which can be extrap¬ 
olated to E> 50 TeV, yielding 


E^ $^[MW] ^ (2.62 ± .20) x 10"® 


E 

100 TeV 


-0.61±.05 


GeV^ fu 


(7) 


per V flavor. The predicted energy flux of the NBR (all flavors) obtained from the ’unatten¬ 
uated GBR’, 


E2 4>^ ^ (1.04 ±0.15) X 10"^ 


E 

100 TeV 


-0.61±.05 


GeV^ fu. 


( 8 ) 


is compared in Fig. 2 to the all-flavors NBR measured with IceGube (Aartsen et al. 2015). 
The separate contributions of the Milky Way (~ 76%) and extragalactic sources (~ 24%) 
to the NBR are shown in Fig. 3. Assuming the all flavors NBR flux to be isotropic, its best 
£t single power-law between 25 TeV and 2.8 PeV by the IceGube collaboration (Aartsen et 
al. 2015), E^ ‘I’nbr ~ (6.69 ± 1.20) x 10“® (E/100 TeV)“°'^°^'°® GeV^ fu is in rough agreement 
with Eq. (8). 


The sky distribution of the NBR measured with IceGube is expected to coincide with 
that of the unattenuated high energy GBR, which is roughly that measured by Fermi-LAT 
at 100 GeV. This distribution that is peaked sharply around the Galactic center is shown 
in Fig. 4. Its peak, however, subtends only a small solid angle: The GBR that was mea¬ 
sured with Fermi-LAT near 100 GeV (Ackermann et al. 2012) suggests that only ~ 4.2% of 
the neutrino events point back towards the Galactic center within latitudes —8° < b < ±8° 
and longitudes —80° < 1 < ±80°, which cover only ^ 0.43% of the full sky. Also plot¬ 
ted are the sky distribution of the GBR at E>1 GeV measured with EGRET aboard the 
Gompton Gamma Ray Observatory and normalized to the flux measured by Fermi-Lat at 
100 GeV. The diffuse gamma-ray emission (flux and sky distribution) from the Galactic 
plane (0° < b < 2°) and from point sources measured at higher energies (e.g., Abramowski 
et al. 2014 and references therein) with buried muon detectors (e.g., GASA-MIA), water 
Gherenkov detectors (e.g., Milagro and HWAG) and atmospheric Gherenkov detectors (e.g., 
H.E.S.S., MAGIG, and VERITAS), are generally consistent within errors with that extrap¬ 
olated from the Fermi-LAT GBR assuming 4*7 (E) oc modified by attenuation in the 

Galactic and extragalactic background light. 


4. Cosmic ray positrons near Earth 

A detailed derivation of the expected flux of high energy cosmic ray positrons near 
Earth produced in hadronic interactions of cosmic ray nucleons in the local ISM, and its 



comparison to the flux measured with high precision by the AMS02 collaboration (Aguilar 
et al. 2015) is presented in Dado & Dar 2015. Here we summarize it for completeness. 

In a steady state, the local flux <l>e+(E) of e’''’s produced in the ISM satisfies 

^lb(E)4,UE)| = J,+ (E) (9) 

where b(E) = —dE/dt is the loss rate of e"*" energy by radiation (rad) and by escape (esc) 
from the Galaxy by diffusion through its turbulent magnetic fields, and 

Je+ (E) Ke+ ain(pp) nism C <hp (10) 

is the local production rate of CR positrons in the ISM whose nucleon density in the solar 
neighborhood is njsm- The solution of Eq. (9) is 

<I>e+(E) ^ Ke+ ain(pp)nismCre<l>p(E)/(/3j - 1) . (11) 

where <hp is given by Eq. (2), njsm ~ 0.9 cm“^ in the solar neighborhood (e.g. Kalberla & 
Dedes 2008), Ke+ ^ 7 x 10“^ for I3j ^ 2.7 — 0.06 = 2.64 and = E/(dE/dt) is the mean life¬ 
time of positrons in the ISM due to their escape from the Galaxy by diffusion (dE/dt E/ Tgsc) 

and radiative energy losses (inverse Gompton scattering of background photons and syn¬ 
chrotron radiation). It satisfies l/rg = 1/resc + 1/rrad, i-e., Tg = Tgsc Trad/(Tesc + TVad) • 

For random Galactic magnetic fields with a Kolmogorov spectrum 

Tgsc ^ 7.5 X 10^^ (E/GeV)-^/^ s (12) 

where the normalization has been adjusted to the value obtained by Lipari (2014) from a 
leaky box model analysis of the flux ratio ^°Be/®Be measured with the Gosmic Ray Isotope 
Spectrometer (GRIS) in the energy range 70-145 MeV/nucleon (Yanasak et al. 2001). 

The radiative life time due to synchrotron emission in the local ISM magnetic held 
with energy density U ~ B^/Stt ~ 0.40eV/cm^ and inverse Gompton scattering of the local 
radiation background (diffuse Galactic light (DGL) with energy density U ~ 0.41eV/cm^, 
far infra red (FIR) light with U ~ 0.40eV/cm^, and cosmic microwave background (GMB) 
with U ~ 0.26eV/cm^), was calculated in the Thomson and Klein Nishina regimes following 
the approximations introduced by Schlickeiser & Ruppel (2010). Other energy loss mecha¬ 
nisms that are important only at energy well below 10 GeV (Goulomb scattering, ionization 
and bremsstrahlung) as well as threshold effects, geomagnetic shielding and solar modula¬ 
tion, were included for completeness through a best fit phenomenological depletion factor 
D(E) = 1 — exp(—(E/V)"), which depends on the time of the measurements but does not 
affect the behavior at E>10 GeV. 
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In Fig. 5, the flux measured with AMS02 (Aguilar et ah 2015) is compared to the 
expected local flux of CR e’''’s as given by Eq.(ll) where Ke+ *hp(E) oc E~^ has been replaced 
by <Fp(E/xe+)/xe+ with the CR proton flux measured with AMS02 (Aguilar et ah 2015). As 
can be seen from Fig. 4, the agreement is quite good. 


5. Conclusions 

We have shown that the observed fluxes of very high energy backgrounds of astronomical 
neutrinos (NBR) and gamma rays (GBR), measured respectively with IceCube and Fermi- 
LAT, satisfy a simple relation, which follows from a common production in high energy 
hadronic collisions of cosmic rays in/near their main Galactic and extragalactic sources. An 
additional stringent test of the common hadronic origin of the high energy NBR and GBR 
is whether the sky distribution of the NBR measured with IceCube is approximately that of 
the unattenuated GBR near 100 GeV. Such a test, however, requires much higher statistics 
than currently available from IceCube. We have also shown that the observed flux of high 
energy cosmic ray positrons near Earth measured with AMS02 is in good agreement with 
that expected from hadronic cosmic ray production in the local ISM. Hence, we conclude that 
at present the observed fluxes, spectra, and sky distributions of the high energy astronomical 
neutrinos, gamma rays and cosmic ray positrons do not provide any compelling evidence for 
decay/annihilation of dark matter particles. 


Appendix A: Photo-pion origin of the NBR and GBR? 

A.l: Upper bound on the NBR ? 

A ’’model-independent upper bound” 

E^2$^ < 2 X 10"® GeV^ fu (13) 

(per neutrino flavor) for the energy flux of extragalactic NBR above 100 TeV was derived 
by Waxman & Bahcall (1999). It was based on the assumption that the NBR is produced 
by the extragalactic cosmic ray protons in proton-photon collisions in GR sources which 
are optically thin. This ’’model independent upper bound”, however is a model dependent 
estimate fraught with incorrect extrapolation (see, e.g., Ahlers, et ah 2005) and incomplete 
consequences of arbitrary choices. When properly corrected, this ’’upper bound” increases 
nearly by two orders of magnitude! and becomes practically useless. 

Waxman & Bahcall (1999) assumed that the flux of ultra high energy (UHE) cosmic 
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rays above the CR ankle observed at Earth is a universal flux of extragalactic CR protons 
with an energy independent injection-rate per comoving volume (Waxman 1995), 

Ecr ^ 10^^ erg Mpc”^ yr"^. (14) 

Ecr 

Assuming that this injection rate is valid all the way down to 100 TeV, and that each 
cosmic ray proton photo-produces at most a single before escaping the source, Waxman 
& Bahcall (1999) obtained a ’’model-independent upper bound” 

< 2 X 10"® GeV^ fu (15) 

(per neutrino flavor) for the energy flux of the extragalactic NBR in the energy range 100 TeV 
to 100 PeV. If the UHE cosmic rays that are observed on Earth are extragalactic in origin, 
they must be protons. This is because of the rather short lifetime of extragalactic UHE CR 
nuclei due to their photo-disintegration in collisions with the extragalactic background light 
(Stecker 1969). A proton-dominated composition is also indicated by the observed break 
near (5 ± 1) x 10^® eV in the spectrum of the UHE cosmic rays observed with the Hires 
(Abbasi et ah 2008), Pierre Auger (e.g., Abraham et al. 2010), and the Telescope Array 
observatories, that coincides with that predicted by Greisen (1966), and Zatsepin & Kuzmin 
(1966) for UHE extragalactic CR protons. However, Waxman & Bahcall (1999) also assumed 
a spectrum dNcn/dE oc E“^ of the UHE extragalactic CR protons, which extends all the way 
down to PeV energies, while the observed flux of the UHE CRs below the GZK cutoff is well 
represented by (Aab et al. 2015) <hp ~ 0.2 (E/GeV)”^ ’^^ fu. This spectral behavior of UHE 
CR protons also well represents, within errors, the measured spectrum of CR protons below 
the all particle CR ankle (Apel et al. 2013) down to ~100 PeV, presumably the CR proton 
ankle (below the proton ankle, diffusion in the Galactic turbulent magnetic fields with a 
Kolmogorov spectrum changes their spectral index to ~ —2.71 — 1/3 ~ —3.04). In a steady 
state, the extragalactic CR spectrum is identical to the CR source spectrum. The observed 
extragalactic flux of CR protons with a power-law index -2.71 (Aab et al. 2015) extrapolated 
down to 100 TeV yields an energy flux 3.5 x 10® times larger than that of an extragalactic 
flux of CR protons with power-law index -2 postulated by Waxman & Bahcall (1999). 

If the production of high energy neutrinos above 100 TeV is dominated by the decay of 
TT®' from p 7 —)■ A'''(1232) nTr®", then the tt^^’s carries on average a fraction E(,./Ma+ ~ 0.22 
of the energy of the incident protons where E(. is the tt energy in the A+ rest frame. Each 
neutrino from the tt’*' decay carries an average energy <Ei,>^E7r/4. Hence, production 
of neutrinos through photo-production of A’'‘(1232)’s by high energy cosmic ray protons in 
optically thin sources satisfies 


^ (1/2) (0.22/4)®-^® Na $p ~ 3.5 x 10^ Na <hp 


( 16 ) 
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per neutrino flavor (assuming complete flavor mixing by neutrino oscillations), where 
Na = o'p 7 1 is the mean number of A+’s produced in collisions of a single UHE CR 

proton with a photon column density before it escapes into the intergalactic space. Since 
the power-law fluxes of CRs and neutrinos from the CR sources have the same power-law 
index, they suffer the same redshift. Hence, for a constant Na per CR proton, all the 
dependence of Eq. (16) on cosmic evolution is through the dependence of the accumulated 
<hp on cosmic evolution. 


But, a high energy CR proton looses ~ 22 % of its energy in the reaction p 7 —)■ A'''(1232) —)■ utt’'' 
followed by n —)■ pe“ Pe- Since the branching ratio for Tr"*" production is ~ 1 / 2 , in order to 
photo produce a single tt’*', each UHE proton must produce on average two A+’s. Hence, 
such protons emerge from their sources with a fraction (0.78)^ ~ 0.61, of their initial 
energy. For a spectral index -2.71, it implies that the in-source flux of such CR protons must 
be larger by a factor (1/0.61)^'’^ ~ 2.34 than that of the emergent flux in order to produce 
a single per observed CR proton. Such a flux produces a neutrino flux larger by a factor 
2.34 than the observed flux of UHE CR protons. Consequently, the maximal neutrino flux 
from ’’thin sources”, which were dehned by Waxman & Bahcall (1999) as sources producing 
a single per CR proton, the energy flux per neutrino flavor produced in-source is 


^ 4.7 < x^ E^ $p(E) ^ 10"® 


E 

100 TeV 


-0.71 


GeV^ fu. 


(17) 


This estimate for ’’maximally thin” sources is larger by nearly two orders of magnitude than 
the Waxman-Bahcall ’’upper bound”, which was based on an assumed power-law index -2 
of the extragalctic CR protons above the CR ankle instead of the observed -2.71 (e.g., Aab 
et al. 2015 and references therein) and was used for extrapolating the energy flux of the 
extragalactic UHE CR protons down to 100 TeV, and which neglected the energy loss of CR 
protons in source in producing ”at least a single tt’*'” (i.e., two A’''’s) per CR proton. 


A.2: Photo-pion origin of the NBR ? 

If the production the NBR and GBR above 30 TeV is dominated by the decay of 
photo-pion produced in CR proton-photon collisions py —)■ A+(1232) —)■ , and —)■ p 7 r°, 

respectively, in/near source (Mannheim 1993), then, for CR protons with in-source spectral 
index [ds = —2.7, the unattenuated NBR and GBR satisfy the relation 

^ (l/ 2 )(m ^+/2 m^o)i-^ ^ 0.16$^ . (18) 

per u flavor. However, if the steepening of the CR spectrum at the CR knee is due to a 
transition from Kolmogorov diffusion to drift motion, while the spectral index of the source 
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spectrum (3^ = —2.3 does not change (Anchordoqui et al. 2014b), then 4),^ 0.214)..^, per v 

flavor. Hence, the extragalactic NBR expected from the observed EGB with Fermi-LAT, 


E2$^[EG] (1.8 ±.30) X 10"^ 


E 

100 TeV 


-0.61±.05 

GeV^ fu. 


(19) 


i.e., much smaller than that measured with IceGube. 


The Galactic contribution of photo-pion production by GR interactions in/near source 
to the NBR below ~ 3 PeV is strongly suppressed by a too high effective energy threshold. 
The effective threshold energy for photo-pion production in collisions of GRs with the typical 
gray body radiation field of a mean photon energy e-^ is 


E. 


'th 


m^ mr 


2e. 


6.5 X 10^6 
eJeV 


eV. 


( 20 ) 


The mean photon energy of stellar light in the ISM and IGM is e..,, ~ 1 eV. Hence, the 
minimal energy of the bulk of neutrinos from Galactic photo-pion production is smaller by 
a factor 20, yielding an energy threshold Ej/ > 3 PeV for Galactic photo-pion neutrinos 
produced by Galactic cosmic rays in/near source and in the ISM. 
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E(MeV) 


Fig. 1.— The flux of the extragalactic gamma-ray background (EGB) as function of gamma- 
ray energy measured with Fermi-LAT (Ackermann et ah 2015) and the best fit exponentially 
cutoff power-law. The straight line represents the unabsorbed power-law EGB. 
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Fig. 2.— Comparison between the energy flux of the NBR (all flavors) above 20 TeV mea¬ 
sured with IceCube (Aartsen et al. 2015) and that expected from the unattenuated GBR 
below 2 TeV as inferred from the GBR and EGB measured with Fermi-LAT (Ackermann 
et al. 2012,2015). Also shown is the best fit single power-law to the unattenuated GBR 
(spectral index -2.36) in the energy range 100 GeV - 2 TeV estimated from the Fermi-LAT 
measurements. 
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Fig. 3.— Comparison between the energy flux of the NBR (per v flavor) above 50 TeV, which 
was inferred from the first 3 years measurements with IceCube (Aartsen et ah 2014), and 
the flux of the NBR (per v flavor) expected from the GBR and EGB, which were measured 
with Fermi-LAT (Ackermann et ah 2012,2015). The separate contributions of extragalactic 
(EG) neutrinos and neutrinos from our Galaxy (MW) to the NBR are also shown. 
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Fig. 4.— The sky distribution of the high energy GBR as function of Galactic latitude, 
observed with Fermi-LAT (Ackermann 2012) at 100 GeV. Also shown is the sky distribution 
observed with EGRET aboard the Gompton Gamma Ray Observatory at E > 1 GeV (Pohl 
et ah 1997) normalized to the Fermi-LAT distribution. The NBR is predicted to have nearly 
the same sky distribution as that of the unattenuated high energy GBR. 
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Energy (GeV) 

Fig. 5.— Comparison between the CR positron flux measured with AMS02 (Aguilar et 
ah 2014) and the flux expected from positron production in the local ISM. 


















